For  Reference 


not  to  be  taken  from  this  room 


Mwnsiww 


University  of  Alberta 
Printing  Department 


Digitized  by  the  Internet  Archive 
in  2019  with  funding  from 
University  of  Alberta  Libraries 


https  ://arch  i  ve .  org/detai  Is/H  rudey  1 967 


THE  UNIVERSITY  OF  ALBERTA 


BENDING  OF  RATE  SENSITIVE 
CIRCULAR  PLATES 


by 


TERRY  MICHAEL  HRUDEY 


A  THESIS 

SUBMITTED  TO  THE  FACULTY  OF  GRADUATE  STUDIES 
IN  PARTIAL  FULFILMENT  OF  THE  REQUIREMENTS  FOR  THE  DEGREE 

OF  MASTER  OF  SCIENCE 


DEPARTMENT  OF  MECHANICAL  ENGINEERING 


EDMONTON,  ALBERTA 
JUNE,  1967 


UNIVERSITY  OF  ALBERTA 


FACULTY  OF  GRADUATE  STUDIES 


The  undersigned  certify  that  they  have  read,  and 
recommend  to  the  Faculty  of  Graduate  Studies  for  acceptance, 
a  thesis  entitled  "BENDING  OF  RATE  SENSITIVE  CIRCULAR  PLATES" 
submitted  by  TERRY  MICHAEL  HRUDEY  in  partial  fulfilment  of 
the  requirements  for  the  degree  of  Master  of  Science. 


ABSTRACT 


A  linearized  theory  based  on  piecewise  flat  level 
surfaces  of  the  complementary  function  is  developed  for 
an  incompressible  inelastic  dissipative  continuum.  Unique¬ 
ness  and  extremum  principles  associated  with  the  constitutive 
equations  derived  from  the  complementary  function  are 
discussed. 

The  theory  is  applied  to  the  problem  of  a  clamped 
circular  visco-plastic  plate  loaded  with  a  uniform  pressure, 
and  the  stationary  creep  of  a  plate  with  the  same  geometry 
and  loading.  Complete  solutions  are  found  for  the  bending 
moments  and  rate  of  deflection  for  both  problems. 
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CHAPTER  I 


INTRODUCTION 


1.1  Preliminary 


The  behavior  of  axially  symmetric  flat  plates  with 
various  edge  conditions  is  of  engineering  importance  and 
has  thus  received  considerable  attention.  Until  recently 
however,  only  the  elastic  problem  had  been  considered. 

Recent  developments  in  continuum  mechanics  have  resulted 
in  solutions  for  this  type  of  plate  problem  for  more  com¬ 
plicated  constitutive  relations. 

Hopkins  and  Prager  9  in  1953  considered  the  problem 
of  a  circular  rigid-perf ectly  plastic  plate  obeying  the 
Tresca  yield  condition.  Closed  form  solutions  were  obtained 
for  the  moments  and  deflection  rates  along  with  the  corres¬ 
ponding  limit  loads  for  both  simply  supported  and  clamped 
edges.  The  solutions  are  valid  only  for  small  deflections 
since  membrane  stresses  are  not  considered.  Also  the  rigid 
plastic  model  is  more  realistic  for  a  sandwich  plate  than 
for  a  solid  plate. 

The  importance  of  viscous  or  rate  dependent  effects 
in  metals  at  elevated  temperatures  has  resulted  in  a 
number  of  attempts  to  take  these  effects  into  account. 

This  has  led  to  the  application  of  various  creep  laws  and 
the  visco-plastic  model  to  plate  problems. 


atrooelv  io  o onelz 
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Venkatraman  and  Hodge 


16 


have  obtained  solutions 


for  creep  of  both  simply  supported  and  clamped  uniformly 

loaded  circular  plates,  A  creep  law  based  upon  a  flow  rule 

associated  with  a  condition  of  maximum  shearing  stress  was 

assumed.  The  solution  assumes  further  that  stationary 

creep  conditions  prevail  and  that  the  creep  rate  can  be 

represented  by  a  function  of  moment  times  a  function  of  time. 

The  assumption  of  stationary  creep  is  equivalent  to  assuming 

that  deformation  occurs  under  a  time  invariant  state  of 
* 

stress  ,  since  changes  in  the  stress  would  result  in  changes 
in  the  elastic  deformation  as  well  as  in  additional  primary 


creep, 


Prager 
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has  suggested  a  linearized  theory  of  visco¬ 
plasticity  which  involves  a  piecewise  linear  yield  condition. 
The  proposed  model  reduces  to  the  Newtonian  viscous  fluid 
when  the  yield  stress  approaches  zero  and  uniqueness  has 
been  proven  for  the  strain  rates.  This  model,  with  the 
Tresca  yield  condition  has  been  applied  to  the  problem  of 
a  simply  supported  circular  plate  under  uniform  pressure 
and  a  complete  solution  has  been  obtained  1  .  It  appears 
that  the  application  of  this  theory  to  the  problem  of  a 
clamped  plate  is  intractable. 


*  In  some  literature  the  term  stationary  creep  is  synonymous 
with  secondary  creep,  that  is  creep  with  constant  strain 
rate.  In  this  thesis  stationary  creep  refers  to  creep 
occurring  under  a  constant  state  of  stress. 


3 


Haddow  3  has  suggested  the  linearized  theory  discussed 
in  this  thesis  and  has  found  a  complete  solution  for  a  visco¬ 
plastic  simply  supported  circular  plate  loaded  with  a  uniform 
pressure „ 

This  thesis  considers  the  problem  of  a  uniformly  loaded 
clamped  circular  plate.  A  linearization  technique  is  dev¬ 
eloped  and  applied  to  two  types  of  isotropic  materials.  The 

first  is  a  rigid  visco-plastic  solid  which  is  a  generalization 

* 

of  the  Bingham  solid.  The  second  is  a  material  obeying  a 
creep  law  which  is  a  generalization  of  the  one  dimensional 
law 

«  -  fe 


.n 


where  and  n  are  material  properties  and  €  and  a  are  the 
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Both  materials 


strain  rate  and  the  stress  respectively 
have  as  special  cases  the  Newtonian  viscous  fluid  and  the 
rigid-perf ectly  plastic  solid.  For  n  =  1  the  creep  law 
above  reduces  to  the  constitutive  equation  in  one  dimension 
for  a  Newtonian  viscous  fluid  and  when  n  approaches  infinity 
the  rigid-perf ectly  plastic  solid  is  approached. 


ie 


The  Bingham  solid  is  a  solid  that  has  the  constitutive 
equation 


12 


=  k  +  2\± 


12 


for  a  state  of  simple  shear  a =  °i2>  w^ere  ^  is  t^ie 
yield  stress^in  pure  shear  and  |i  is  a  coefficient  of 
viscosity 
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1.2  The  Work  and  Complementary  Functions 

Hill  5  has  attempted  to  establish  a  unified  theory 
for  several  branches  of  the  subject  of  continuum  mechanics. 
The  basis  of  this  unified  theory  depends  on  the  existence 
of  two  convex  functions,  the  work  function  E  and  the  comple¬ 
mentary  function  Ec . 

Let  o. .  and  e . .  be  the  components  of  the  stress  and 
ID  ID 

strain  rate  tensors  respectively,  referred  to  the  rectangular 

axes  ox. .  If  a .  .  is  a  sinqle  valued  function  of  e .  .  then 
1  ID  ID 

for  a  purely  dissipative  medium  the  stress  tensor  can  be 
written  as 


o .  .  = 

ID 


dE 

de  .  . 
ID 


(1.1) 


where  for  the  purpose  of  differentiation  all  nine  components 
of  e^j  are  considered  to  be  independent.  The  function  E(e^j) 
is  the  work  function  and  is  defined  as 


E(eij>  = 


a .  .  de .  . 
ID  id 


(1.2) 


The  condition  that  E  be  a  convex  function  is 
*  .  -  *  dE 


<eij  >  -  E(eij>  >  (eij  -  eij>  5eT 


ID 


for  any  two  strain  rate  states  e^^  and  e^j.  Considering 

* 

a  surface  of  constant  E  in  principal  strain  rate  space  , 


(1.3) 


The  concept  of  principal  strain  rate  space  or  principal 
stress  space  is  useful  only  for  isotropic  media.  For 
anisotropic  media  9  dimensional  spaces  must  be  considered. 


■ 
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where  the  co-ordinate  axes  are  the  three  principal  strain 
rates,  the  condition  (1«3)  requires  that  the  surface  is  no¬ 
where  convex  to  the  origin 0  If  the  strict  inequality  in 
(1.3)  is  always  valid  when  e^j  is  not  equal  to  e^j  then 
the  surface  is  said  to  be  strictly  convex,  meaning  that 
a  surface  of  constant  E  is  everywhere  concave  to  the  origin 
and  contains  no  flat  surfaces,, 


The  complementary  function  Ec  is  defined  in  a  similar 

manner.  If  e. .  is  a  single  valued  function  of  o. .  then  for 

1  j  ^  1  j 

a  purely  dissipative  medium  the  components  of  the  strain 
rate  may  be  written  as 

dE 


ID 


13 


where  is  called  the  complementary  function  and  is 


defined  as 


E  (o.  . 

c  13 


e .  .  do .  . 
ID  ID 


(1.4) 


(1.5] 


The  condition  for  E  (o.  .)  to  be  convex  is  analagous 

c  ij 

to  (1.3)  and  is 


*  *  dE 

E  (o.  .  )  -  E  (o.  .)  >  (a.  .  -  a.  .) 

c  13  c'  13  13  13  00 


ID 


for  any  two  stress  states  o .  .  and  a 

■*  1313 


* 


From  the  definitions  of  E  and  Ec  it  follows  that  the 

existence  of  both  functions  implies  a  one-one  relationship 

between  o. .  and  e. .  if  both  E  and  E  are  strictly  convex. 
13  13  c 


(1.6) 


' 
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The  sum  of  the  work  and  complementary  functions  is 


E  +  E 


a .  .  de .  .  +  /  e  .  .  da .  . 

ID  ID  J  ID  ID 


e  .  .  a  .  . 
ID  id 


(1.7) 


and  is  the  rate  of  energy  dissipation  per  unit  volume. 

If  all  nine  components  of  e^j  and  o ^  are  considered 
as  independent  variables,  equations  (1.1),  (1.4),  and  (1.7) 


are  seen  to  be  a  Legendre  dual  transformation 


10 


If  the  strain  rate  components  e^^  are  taken  as  co¬ 
ordinates  in  a  nine  dimensional  space,  equation  (1.1)  indi¬ 
cates  that  the  are  components  of  a  vector  normal  to  a 

surface  of  constant  E,  the  components  of  a^^  being  multi¬ 
plied  by  a  scalar  constant  to  give  them  the  dimensions  of 
strain.  Similarly  from  equation  (1.4)  the  e^  are  seen  to 
be  the  components  of  a  vector  normal  to  a  surface  of  con¬ 
stant  Ec  in  nine  dimensional  stress  space. 

If  the  material  is  incompressible,  the  first  invariant 

of  the  strain  rate  tensor,  ev ,  ,  is  zero.  Consequently 
c)E_  KK 

-  is  zero  and  E  is  independent  of  the  hydrostatic  com- 


kk 


ponent  of  stress  a. 


.  i  .  Thus  the  stress  tensor  a.  .  can  be 
kk  i  j 


replaced  by  its  deviator  s^j,  where 


s  .  . 
ID 


a .  . 

iD 


a  5  .  . 

3  kk  ID 


(1.8) 


An  example  of  an  incompressible  material  for  which  the 
work  and  complementary  functions  exist  is  the  classical 
Newtonian  viscous  fluid.  The  two  convex  functions  for  this 


material  are 


. 


7 


E 


H  e 


ID 


J 

2m 


(1.9) 

(1.10) 


where  I  -  2e . .  e. .  and  J  =  ~  s.  .  s. .  are  invariants  of  the 
ID  ID  2  13  13 

strain  rate  and  stress  deviator  tensors  respectively,  and  m 
is  the  viscosity.  The  constitutive  relation  which  may  be 
derived  from  either  (1.9)  or  (1.10)  is 


(1.11) 


Equation  (1.10)  is  a  special  case  of  a  more  general  complement¬ 
ary  function  of  the  form 


Ec  =  f(J) 


(1.12) 


If  the  principal  stresses  o^,  a^,  and  are  taken  as 
rectangular  co-ordinate  axes,  a  surface  of  constant  Ec  given 
by  (1.12)  is  a  circular  cylinder  whose  axis  makes  equal  angles 
with  each  of  the  co-ordinate  axes. 

A  further  example  of  a  material  with  a  complementary 
function  of  the  form  given  by  (1.12)  is  the  generalization  of 

*  r  " 


the  Bingham  solid  suggested  by  Prager 


and  Prager 


14 


and  Hohenemser 


for  which 


-  (Ji- 


K>  2/2u 


(1.13) 


*  Often  in  the  literature  this  solid  is  referred  to  as  the 
Bingham  solid.  Actually  it  is  a  generalization  of  the  one 
dimensional  model  suggested  by  Bingham.  Henceforth  it  will 
be  referred  to  in  this  thesis  as  the  generalized  Bingham  solid. 
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and  E  -  \±  j;  +  K  JT~  I/O  (1.14) 

where  p,  is  the  viscosity  and  K  is  the  yield  stress  in  pure 
shear  for  the  Mises  yield  condition .  The  angle  brackets 
have  the  following  significance. 

<(f)>  =  0  for  F  <  0 

=  F  for  F  0 


The  constitutive  equations  derivable  from  the  work  and 
complementary  functions  by  use  of  (1.1)  and  (1.4)  are 


and 


1 D 


13 


2(\i  + 


2n  (1 


K 


. — )  e  .  . 

JT 

K  \ 

"VF/ 


I/O 


(1.15) 

(1.16) 


From  (1.15)  it  follows  that  the  stress  in  the  generalized 

Bingham  solid  subjected  to  a  strain  rate  e. .  is  the  sum  of  the 

1 3 

stresses  in  a  Newtonian  viscous  fluid  and  a  Mises  rigid  plastic 
solid  each  subjected  to  the  same  strain  rate. 

The  work  function  (1.14)  is  equal  to  one  half  the  rate 
of  viscous  energy  dissipation  per  unit  volume  plus  the  total 
rate  of  plastic  energy  dissipation  per  unit  volume.  The 
complementary  function  (1.13)  is  equal  to  one  half  the  rate 
of  viscous  energy  dissipation  per  unit  volume. 

The  complementary  and  work  functions  for  both  the 
Newtonian  viscous  fluid  and  the  Bingham  solid  are  both  strictly 
convex  and  thus  the  uniqueness  theorems  due  to  Hill  are  valid. 
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1.3  Linearization 


The  basic  hypothesis  of  the  classical  inviscid  theory 


of  plasticity 


is  the  existence  of  a  yield  condition 


f(cKj)  =  c  and  the  equivalence  of  the  yield  function  f  with 
the  plastic  potential.  The  plastic  strain  rates  e?j  for  a 
rigid-perf ectly  plastic  material  are  thus  given  by 


,  P  = 
ij 


A 


df 


13 


where  A  is  a  non-negative  parameter  and 


(1.17) 


A  =  0  if  f  <  c  or  if  f  =  0  and  f  <C  0 

A  ^  0  if  f  =  c  and  f  =  0 


where  f  is  the  material  derivative  of  f. 

In  principal  stress  space  the  vector  with  components 
(e^,  e 2 ,  e^)  is  normal  to  the  yield  surface  f(cKj)  “  c •  For 
the  Mises  yield  condition 


f  (0.  .) 

13 


(1.18) 


9 

and  c  =  k"  where  k  is  the  yield  stress  in  pure  shear.  The 
yield  surface  is  a  circular  cylinder  whose  axis  passes  through 
the  origin  and  makes  equal  angles  with  each  of  the  co-ordinate 
axes . 


A  similarity  in  form  can  thus  be  seen  between  the  flow 
rule  for  a  Mises  rigid  plastic  solid  and  a  material  with  a 
complementary  function  given  by  (1.12). 


. 
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The  Tresca  yield  condition  is  given  by 

la  -  a  .  |  =  2k 

I  max  mm  I 


(1.19) 


where  a _ and  a  .  are  the  maximum  and  minimum  principal 

max  mm 

stresses  respectively  and  k  is  the  yield  stress  in  pure 
shear.  The  yield  surface  for  a  Tresca  solid  is  a  regular 
hexagon  in  principal  stress  space.  It  may  further  be  re¬ 
garded  as  a  piecewise  linear  approximation  for  the  Mises 
yield  surface. 

This  suggests  the  possibility  of  a  similar  piecewise 

approximation  for  a  complementary  function  of  the  form  (1„12)0 

* 

Such  an  approximation  could  be  given  in  general  as 


Ec  =  g(G) 


(1,20) 


where 


o  -  a  -  a  . 

max  mm 


(1.21) 


By  a  choice  of  the  proper  function  g(a) ,  the  surface  of 
constant  Ec  given  by  (1,20)  could  be  made  to  fit  inside  the 
surface  of  constant  Ec  given  by  (1,12),  such  that  the  two 
surf aces co-incide  only  at  the  corners  of  the  hexagonal 
prism  (1,21) , 

This  approach  when  applied  to  the  Newtonian  viscous 


*  Although  this  is  an  approximation  in  the  sense  that  a  cylin¬ 
drical  surface  of  constant  E  is  replaced  by  a  hexagonal  prism 
in  principal  stress  space,  tfiere  is  experimental  evidence  17 
that  the  constitutive  equations  derived  from  the  latter  [_ 
form  of  the  Ec  surface  represent  more  closely  the  creep  behavior 
of  some  materials.  This  is  analagous  to  the  fact  that  the  yield 
condition  for  most  real  metals  lies  somewhere  between  the 
Mises  and  the  Tresca  yield  conditions. 


11 


fluid  yields  the  following  form  of  the  complementary  function. 


/  \  2 
(a  -  a  .  ) 
max  mm 


(1.22) 


The  intersection  of  the  two  surfaces  of  constant  E 

c 

given  by  (1.10)  and  (1.22)  with  the  plane  ~  constant, 
in  principal  stress  space,  are  shown  in  Fig.  1.1.  The  posi¬ 
tion  of  the  point  (o^,  is  determined  by  the  hydrostatic 

component  of  stress. 


FIGURE  1.1 

SURFACES  OF  CONSTANT  E  FOR  VISCOUS  FLUID 

c 

If  the  stress  point  lies  on  a  flat  of  the  surface  of 
constant  Ec ,  the  strain  rate  vector  corresponding  to  that 
state  of  stress  is  normal  to  the  flat.  Thus,  if  the  three 
principal  stresses  are  distinct,  meaning  that  the  stress 


12 


point  is  on  a  flat,  the  principal  components  of  the  strain 
rate  obtained  from  equation  (1.4)  are 


max 


=  (a 


max 


o  •  )  /3|i 
mm 


'min  ^ 


int 


=  0 


(1 


where  e  ,  e.  , ,  and  e  .  are  the  maximum,  intermediate, 
and  minimum  principal  components.  Since  this  is  true  for 
any  point  on  the  flat  the  components  of  the  stress  tensor 
are  not  single  valued  functions  of  the  components  of  the 
strain  rate  for  a  material  with  this  form  of  the  complementary 
function . 

If  any  two  of  the  principal  components  of  stress  are 

equal  the  stress  point  lies  on  a  corner  of  the  surface  of 

constant  E  .  The  strain  rate  vector  which  is  normal  to 
c 

the  surface  of  Ec  is  thus  not  uniquely  determined.  The 
strain  rate  vector  corresponding  to  such  a  stress  state, 
can  lie  anywhere  within  the  fan  enclosed  by  the  normals  to 
the  adjacent  flats  of  the  Ec  surface.  Although  the  direction 
of  the  strain  rate  vector  is  not  unique  there  is  a  restriction 
on  the  magnitude  of  the  vector.  The  three  principal  com¬ 
ponents  of  the  strain  rate  vector  for  any  vector  lying  in 
the  fan  enclosed  by  the  normals  to  the  two  adjacent  flats 

of  the  E  surface,  must  satisfy  the  condition  that  the  rate 
c 

of  energy  dissipation  per  unit  volume  is  constant  for  any 
point  on  the  Ec  surface.  This  follows  from  the  fact  that 
for  this  material,  Ec  is  one  half  of  the  total  rate  of  energy 


23) 


nxsi^a  arf* 
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dissipation  per  unit  volume.  This  condition  may  be  stated 
as 

+0262+  CJ3  =  constant  .  (1.24) 

Referring  to  Fig.  1.1,  the  strain  rate  vector  for  a 
point  on  the  flat  AB  has  components 

(°1  -  °3> 

el  =  - 3ir^-  =  -e3  *  e2  =  0  (1-25) 

and  the  strain  rate  vector  for  a  stress  point  on  the  flat 
CB  has  components 


The  strain  rate  vector  for  point  B  can  thus  be  found 
using  (1.24),  (1.25),  and  (1.26).  It  has  the  components 


t  (c?!  -  g3)  (1-t)  (o2  -  o 3) 

el  =  3l T  e2  =  3'ul  *  e3  =  "el  ”  e2  (!U27) 

where  t  is  a  parameter  which  can  take  any  value  in  the  closed 
interval  0^  t  ^1.  Similar  expressions  may  be  found  for 
each  point  on  the  Ec  surface. 

It  should  be  noted  that  the  strain  rates  resulting  from 
this  approximation  of  Ec  are  independent  of  the  intermediate 
principal  stress.  This  is  also  true  for  the  Tresca  yield 
condition  and  its  associated  flow  rule. 

This  linearization  technique  can  be  applied  to  materials 
with  more  complicated  constitutive  relations.  Consider  now 


14 


a  material  in  which  the  stress  is  given  by  the  sum  of  the 
stresses  in  a  Tresca  rigid-perf ectly  plastic  solid  and  a 
viscous  fluid  with  a  complementary  function  given  by  (1.22), 
when  all  three  bodies  are  subjected  to  the  same  strain  rate. 
This  model  is  a  generalization  of  the  Bingham  solid.  The 
complementary  function  for  this  material  is  found  to  be 


<c 


max 


a 


mm 

6\i 


(1.28) 


where  |jt  is  a  coefficient  of  viscosity  and  k  is  the  yield 
stress  in  pure  shear  for  the  Tresca  yield  condition.  Again 
the  surface  of  constant  Ec  is  a  regular  hexagonal  prism  in 
principal  stress  space. 

In  a  manner  similar  to  that  used  previously  for  the 
viscous  fluid,  the  strain  rates  may  be  obtained.  For  a 
stress  point  on  a  flat  of  the  Ec  surface  the  principal  strain 
rate  components  are 


e 

max 


/a 

\  i 


max 


o 


mm 
3  d 


-e  .  ,  e.  .  =  0  . 

mm’  int 


(1.29) 


As  before,  the  strain  rate  vector  for  a  stress  point 

on  a  corner  of  the  E  surface  must  lie  in  the  fan  enclosed 

c 

by  the  normals  to  the  two  adjacent  flats,  with  the  restriction 
(1.24)  on  the  magnitude  of  the  strain  rate  vector. 

Fig.  1.2  shows  the  intersection  of  the  surfaces  of 
constant  Ec  given  by  (1.13)  and  (1.28),  and  the  plane 
o ^  ~  constant,  along  with  the  intersection  of  the  Tresca 


’ 
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and  Mises  yield  surfaces  with  this  plane. 


1.28) 


FIGURE  1.2 

SURFACES  OF  CONSTANT  E  FOR  VISCO-PLASTIC  MATERIAL 

c 

Table  I  gives  the  strain  rates  corresponding  to  any 
stress  state  for  this  material  and  thus  constitutes  a 
complete  flow  rule. 
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Stress 

Regime 

BOA 

OB 


BOG 

OC 


COD 

OD 


TABLE  I 

FLOW  RULE  FROM  EQUATION  (1.29) 
(Referred  to  Fig.  1.2) 


Stress 

Components 

°1  >  °2  >°3 


Strain  Rate  Components 


(1-t)<°2-°3 


3  \± 


j 


°2  >  °1  >  °3 


-e 


2 


°2  >  °3  >  °1 


e3  el» 

0  <  s  <1 


=  0 


°3>  °2>  °1 


3|± 


DOE 
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Stress  Stress 

Regime  Components 


OE  °3  >  °2  =  01 


Strain  Rate  Components 


(m-1)  <o3”01-2]^> 


:1 


3  |JL 


e2  = 


■m  <^-c 


3  2 


3u 


e3  e2 


■e^  ,  0  m  ^  1 


EOF  °3  ^  °l'>  °2 


el  ^ >  e2  ”e3 ’  e3 


<a3-o2-2k> 


3  _2_ 

3n 


OF  0~,  = 


=  °1> 


a. 


e,  = 


e-,  = 


g<qi-g2-2K> 

3|-l 

(g-l)<a  -o,-2k> 


e2  e3“el 


3il 


,  0  <g  <1 


FOA  01  >  03  >  c2 


OA  a1  =  o. 


3  2 


e,  - 


<°i 


°2‘2k> 


3m- 


-  "e2“e3  ?  e2  ~ 
_  -f<°l-g3-2k) 

3u 


e2  ”ei>  e3  ® 

(f-l)  (o1-a2-2k) 

* 

,  0<f <1 


CHAPTER  II 


UNIQUENESS  AND  EXTREMUM  PRINCIPLES 


2  „  1  Uniqueness 


The  classical  boundary  value  problem  for  a  dissipative 
inelastic  medium  is  posed  as  follows.  Determine  the  distri¬ 
bution  of  stress  and  velocity  if  the  surface  traction  is 

prescribed  on  part  S^  of  the  surface  of  a  body  and  the  vel- 

* 

ocity  v^  is  prescribed  on  the  part  S^  of  the  surface.  The 
uniqueness  of  the  solution  for  such  a  problem  depends  on 
the  nature  of  the  constitutive  relationship  for  the  material. 

The  constitutive  equations  of  the  continua  which  are 
considered  in  this  thesis,  except  for  the  rigid-perf ectly 
plastic  solid,  are  obtained  from  the  complementary  function 
E  .  The  two  types  of  continua  considered  have  as  special 
cases  the  Tresca  rigid-perf ectly  plastic  solid  and  the  viscous 
fluid  with  a  complementary  function  given  by  (1.22). 

The  uniqueness  theorem  for  a  rigid  plastic  solid  has 


been  given  by  Hill 


If  (a .  . , 

i  j  s 


two 


vi)  and  (oij,  v±)  are 
solutions  for  the  rigid  plastic  boundary  value  problem,  then 
by  the  virtual  work,  equation 


*  Surfaces  and  Sv  may  overlap.  Then  the  force  (velocity) 
normal  to  the  surface  and  the  velocity  (force)  tangent  to 
the  surface  are  prescribed  in  the  overlapping  region. 

**  This  uniqueness  theorem  and  the  corresponding  extremum 

principles  are  valid  for  a  rigid  plastic  workhardening  solid 
provided  the  current  distribution  of  hardening  is  known. 
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* 


(F.-F.) (v.-v.)dS  =  7  (a 


* 


.  .-a .  . 

13  13 


13  13 


r  * 

.  _ 

j  (k-x  ) 

V 

dS 


v 


D 


(k-T) 


* 


V 


dS 


* 

D 


D 


(2.1) 


D 


where 


and 


°k 


V 


are  the  velocity  discontinuities  on  S 


D 


and  respectively  and  t  and  T  are  the  shear  stresses  cor- 

responding  to  o. .  and  a. .  acting  on  S  and  S  in  the  directions 

13  X  J  L)  D 


k 


V 


and 


v 


respectively, 


of  the  velocity  discontinuities 

To  show  uniqueness  it  is  necessary  to  introduce  the 
maximum  work  inequality  for  a  rigid  plastic  material, 


13 


* 


13' 


e  ■  •  ^  0 
13  ^ 


2 . 2' 


where  e^  is  the  strain  rate  corresponding  to  the  stress 

* 

in  a  plastically  deforming  body  and  ck  ^  is  any  stress  point 

lying  on  or  within  the  current  yield  surface.  If  the  yield 

surface  is  strictly  convex  then  the  equality  in  (2,2)  can 

* 

hold  only  if  and  are  equal.  If  the  yield  surface 

contains  flats  then  the  equality  in  equation  (2„2)  is  true 
only  if  the  two  stress  points  lie  on  the  same  flat  of  the 
yield  surface. 

ck  ck 

Since  F.  =  F .  on  S,  and  v.  =  v.  on  S  ,  the  left  hand  side 
lif  1  1  v* 

of  (2.1)  must  equal  zero.  From  the  maximum  work  inequality 
and  the  inequalities 

k  >  T,  k  >  T* 


where  k  is  the  yield  stress  in  pure  shear,  it  follows  that 
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each  term  on  the  right  hand  side  of  (2.1)  is  positive  and 
therefore  each  integrand  is  zero.  That  is 


* 

a.  .) 

ID 


0  and 


(2.3) 


Thus  if  the  yield  surface  is  strictly  convex,  equations 
(2.2)  and  (2.3)  indicate  that  the  stress  field  is  unique 
except  perhaps  in  the  common  non-deforming  region  of  the  two 

solutions.  If  the  yield  surface  has  flats  then  equations 

* 

(2.2)  and  (2.3)  require  that  a.  .  and  a.  .  be  on  the  same  flat 

1 D  1 D 

everywhere  except  perhaps  in  the  common  non-def orming  region. 

It  follows  that  the  stress  solution  in  the  deforming  region 
is  not  necessarily  unique.  The  velocity  solution  is  not  unique 
regardless  of  whether  the  yield  surface  has  flats  or  not. 

In  the  axially  symmetric  Tresca  rigid-perf ect ly  plastic 
plate  problem,  which  is  a  special  case  of  the  two  problems 
in  Chapters  III  and  IV,  an  a  priori  assumption  is  made  as  to 
the  stress  regimes  in  which  the  stress  solution  lies  and  a 


solution  is  obtained 
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The  uniqueness  theorem  for  a 
Tresca  rigid  plastic  material  indicates  that  all  stress  solu¬ 
tions  must  lie  on  the  same  stress  profile.  It  is  found  that 
for  any  particular  stress  profile,  only  one  stress  solution 
at  most  can  be  found  which  satisfies  the  equilibrium  equation, 
the  flow  rule,  and  the  boundary  conditions.  Thus  the  stress 
solution  which  has  been  found  elsewhere  for  the  circular 
clamped  Tresca  rigid-perf ect ly  plastic  plate  problem  is 
unique 


# 

' 
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It  remains  now  to  investigate  the  boundary  value 
problem  for  a  material  with  a  constitutive  relation  of  the 
general  form 

dE 


e 


(2.4) 


If  the  stress  in  the  material  is  dependent  on  the  time 
rate  of  deformation  then  velocity  discontinuities  cannot 
occur  since  infinite  stresses  would  be  required  to  produce 
such  discontinuities,,  Therefore  the  virtual  work  equation 
becomes 


* 


13' 


dV  =  0 


(2.5) 


where  (o. e. .)  and  (a. .  ,  e. 

13  13  13s  13' 

boundary  value  problem . 


are  two  solutions  to  the 


Substitution  of  the  constitutive  equation  ( 2 » 4)  into 
( 2  o  5)  gives 


v 


* 

o  .  .  -  o  .  . 

1  j  1.  j 


IE  (a . 
c  it 

13 


dV 


0 


(2.6) 


Equation  (2.6)  is  satisfied  if 
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13 


*  dE  (a..) 

a*  . )  —  S— 13-  = 


13 


1  o 


(a 


13 


13 


a 


dE  (o.  .) 

C  IT 


13 


do 
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From  the  condition  of  convexity  for  the  Ec  surface  the 
following  inequality  is  obtained 


G 


2JL 


o .  . 

13 


*  * 

E  (  0  .  .  )  A  ( G  .  .  -  0  .  . 
c '  it  ^  1 3  i  j 


dE  (o.  .) 
c  13' 
* 

do  .  . 

13 


(2.8) 
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If  the  surface  of  constant  is  strictly  convex  then 
the  equality  in  (2.8)  holds  only  if  cs\  ^  and  CK  ^  are  equal. 

It  follows  therefore  from  (2.6)  and  (2.7)  that  the  stress 
solution  is  unique.  Furthermore,  since  Ec  is  strictly  convex 
there  is  a  one-one  relationship  between  stress  and  strain  rate 
and  thus  the  strain  rate  is  also  unique. 

If  E^  is  not  strictly  convex  then  the  equation  (2.7) 

* 

and  the  equality  in  (2.8)  are  satisfied  if  ck  ^  and  o^.  lie 
on  the  same  flat  of  the  surface  of  constant  Ep,  This  implies 
that  the  stress  profile  is  unique. 

The  materials  considered  in  the  two  plate  problems  in 
this  thesis  have  complementary  functions  whose  surfaces  are 
hexagonal  prisms  in  principal  stress  space.  For  these  prob¬ 
lems  the  stress  profile  is  thus  unique.  As  for  the  Tresca 
rigid-perf ectly  plastic  plate  problems  it  is  found  that  for 
any  stress  profile  only  one  solution  can  be  found  which 
satisfies  the  equilibrium  equation ,  flow  rule,  and  the  bound¬ 
ary  conditions.  It  follows  that  the  complete  solution  for  a 
plate  problem  with  a  constitutive  relation  defined  in  terms 
of  the  complementary  function  is  unique. 


*  The  surfaces  of  constant  Ec  considered  in  this  discussion 
are  assumed  not  to  be  convex  to  the  origin  at  any  point . 

Thus  the  distinction  between  a  strictly  convex,  and  a  convex 
E  surface  is  that  the  latter  may  have  flats. 
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2 „ 2  Extremum  Principles 

Two  extremum  principles  may  be  derived  for  a  dissipa¬ 
tive  material  for  which  the  complementary  and  work  functions 
exist.  They  have  been  presented  but  not  proved  by  Hill 

If  (o. v.)  are  the  solution  stress  and  velocity  fields 

ij J  i 

* 

for  this  classical  boundary  value  problem  and  ck  ^  is  any 
statically  admissible  stress  field  then  it  follows  from  the 
virtual  work  equation  that 
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F  .  )  v  .  dS 

i  i 


o  .  .  -  a .  . )  e  .  .  dV . 

13  13  ID 
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But  on  and  from  the  convexity  of  the  com¬ 


plementary  function 


2,9) 
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( a  .  .  -  a  .  . )  e .  .  ^  E  ( a  .  . )  -  E  (a.. 

v  ID  id  id^  c'  id  c  id 


Thus  equation  (2.9)  becomes 
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F  .  v .  dS 
1  1 


E  (0. .)  dV 

c  id 


v.  dS 
1  1 
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Ec(0ij)  dV 
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(2.10) 


which  is  the  first  extremum  principle. 

If  (o.  v.)  are  the  solution  stress  and  velocity  fields 

1  ”1  1 


and  v^  is  any  kinematically  admissible  velocity  field  for 
this  boundary  value  problem  then  from  the  virtual  work 
equation 

f  Fi(vi  ”  vi)  dS  +  f  (vi  -  v^)  dV  =  T  "  eij)dV  (2,11) 


v 


v 


where  X.  is  the  body  force  per  unit  volume, 
x 


. 


s 1  -30  r»v  i>as  noi^oloa  ^h*j  9735  (  ^.v  ,  ._  0)  il 
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* 

On  =  v^s  and  from  the  convexity  of  the  work 

function 


o 


il 


<E 


11 


E  (e 


'll' 


o 


Thus  equation  (2,11)  becomes  the  second  extremim 
principle 


E  ( e 


il 


dV 


(2.12) 


The  right  hand  sides  of  both  (2.10)  and  (2.12)  are 
necessarily  unique. 


CHAPTER  III 


RIGID  VISCO-PLASTIC  PLATE  PROBLEM 

Let  r,  ©,  and  z  be  cylindrical  co-ordinates  with  z 
taken  positive  downwards ,  Consider  a  plate  bounded  by  the 
surfaces  z  =  +  h/2  and  r  =  a, 

The  plate  is  clamped  along  the  edge  r  =  a  and  is  loaded 
with  a  uniform  pressure  p  acting  in  the  positive  z  direction. 
The  rigid  plastic  or  rigid  visco-plastic  model  is  realistic 
for  a  sandwich  type  construction  consisting  of  two  thin  outer 
layers  which  carry  stresses  in  the  plane  of  the  plate  only, 
and  an  inner  layer  which  is  capable  of  carrying  vertical  shear 
only.  When  bending  occurs  these  outer  layers  are  in  a  bi¬ 
axial  state  of  stress.  Since  membrane  stresses,  resulting 
from  extension  of  the  middle  surface  of  the  plate,  are 
assumed  small  in  comparison  to  the  bending  stresses  and 
since  the  stress  is  essentially  constant  across  the  thickness 
of  these  outer  layers,  the  bending  moments  per  unit  length 

M  and  Mn  are  proportional  to  the  stresses  in  the  outer 
r  y 

layers.  Thus  the  generalized  stresses  M^_  and  are  used 
in  the  following  discussion  rather  than  and  . 

Let  w  be  the  velocity  of  any  point  in  the  middle  surface 
of  the  plate,  taken  positive  in  the  positive  z  direction. 

The  usual  assumptions  in  small  deflection  plate  theory  give 
the  curvature  rates  in  the  radial  and  circumferential 


directions  as 


. 

' 


and 
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K. 

r 


_  1  dw 

K0  r  dr  • 


(3.1) 


The  curvature  rates  and  the  moments  are  taken  to  be 
positive  if  they  correspond  to  compression  in  the  z  =  -h/2 
surface.  These  curvature  rates  may  be  regarded  as  general¬ 
ized  strain  rates. 

Consider  now  an  element  of  the  plate  bounded  by  the 
circles  with  radii  r  and  r+dr  and  two  radial  lines  contain¬ 
ing  an  elemental  angle  d@ .  Let  Q  be  the  vertical  shear 
stress  per  unit  length  acting  on  the  surfaces  r  =  constant 
as  shown  in  Fig.  3.1. 


dM 
_ r 

dr 


dr 


FIGURE  3.1 


PLATE  ELEMENT 
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Summing  forces  in  the  z  direction  results  in  the 
following  equilibrium  equation 

+  prdrd©  -  (Q  +  dr)  (r  +  dr)  d©  +  Qrd©  =  0  .  (3.2) 


Neglecting  products  of  differentials  of  degree  greater 
than  two  yields  the  following  equation , 


d(rQ) 

dr 


pr  =  0 


(3,3) 


Equilibrium  of  the  moments  acting  about  a  radial  line 
through  the  element  is  automatically  satisfied  due  to  the 
circular  symmetry  of  the  problem. 

Summing  moments  about  the  circumferential  direction 
gives  the  following, 

dM 

-ML  drd©  +  (M  +  — ~  dr)  (r  +  dr)  d©  -  M  rd©  +  Qrdrd©  =  0.  (3.4) 

©  r  dr  r 


Neglecting  higher  order  small  quantities  this  reduces 


to 


d(rM  ) 
r 

dr 


-  Mq  +  Qr  =  0 


(3.5) 


Integrating  equation  (3.3)  from  0  to  r  gives 


=  ££ 


3.6) 


Substitution  of  this  into  equation  (3.5)  yields  the 
complete  equilibrium  condition  for  the  plate  in  terms  of 


the  generalized  stresses  and  M 


© 


d(rM  ) 
r 

dr 


M  +  2JL 
u©  2 


=  0 


(3.7) 


28 


The  constitutive  relation  which  is  to  be  applied  to 
this  plate  problem  is  defined  in  terms  of  a  complementary 
function  of  the  form  (1.28)  and  may  be  regarded  as  an  approxi¬ 
mation  to  the  generalized  Bingham  solid.  The  complementary 

function  may  be  expressed  in  terms  of  the  generalized  stresses 

2 


M  and  as 
r  © 


E 


<M  -  Mo) 

2  A 


where  M  =  max  (  |  M^_  |  ,  |Mq|  ,  |m^  -  Mq  |  )  and  A  is  a  coefficient  of 

viscosity  whose  relation  to  (i  depends  on  the  plate  geometry. 

The  complementary  function  given  by  equation  (3.8)  is 

equal  to  half  the  rate  of  viscous  energy  dissipation  per 

unit  surface  area  of  the  plate.  If  the  middle  surfaces  of 

the  two  outer  layers  of  thickness  t  are  separated  by  a  dis- 

2  2 

tance  h  then  A  =  3|it  h  .  Similarly  the  yield  moment  is 
related  to  ks  the  yield  stress  in  pure  shear  according  to 
Mq  =  2kth. 

The  curvature  rates  can  be  found  by  partial  differen¬ 
tiation  of  the  complementary  function  with  respect  to  the 


moments  M  and  M 
r  u 


This  results  in  the  flow  rule 


K. 


bE 

c 

cTm” 


dE 


and 


K. 


© 


cM 


(3.9) 


© 


Referring  to  Fig.  3.2  which  shows  a  surface  of  constant 
Ec  given  by  equation  (3.8) ,  the  flow  rule  in  Table  II  is 
obtained  from  equations  (3.8)  and  (3.9). 


. 
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FIGURE  3.2 

SURFACE  OF  CONSTANT  E 

c 

TABLE  II 

FLOW  RULE  FROM  EQUATIONS  (3.8)  AND  (3.9) 


Stress 

Regime 

Complementary 

Function 

Flow  Rule 

COB 

<Me-Mo>2 

'  2  A 

A/c  = 

r 

0 

Ako 

<M0  -  Mq> 

OC 

<Me-Mo>2 

2  A 

>* 

?s 

II 

-s^^-Mq')  0  <  s  ^  1 

Ak9  = 

(l-s)  (Me-M^+s<Me-Mr-M0) 

• 
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Stress  Complementary 
Regime  Function 


COD 


/m  -M  -M 

V  Q  r  0/ 

2A 


Flow  Rule 
X,r  -  -  <Me-Mr-M0> 

X(cQ  -  <M0-Mr-Mo^ 

A(cr  =  -t  <-Mr-Mo>  +(t-l)<M0-Mr-M^> 
XK@  -  (1-t)  <Me-Mr-M^>  0  <  t  ^  1 

XKr  “  "<-Mr-M0> 


0  <  q  <1 


A/c 

OF 

<-M@-Mo>2 

2  A 

A/c 

A/c 

FOA 

<Mr-M0-Mo\2 

2  A 

A/c 

A/c 

OA 

<Mr-M0>2 

2  A 

A/c 

9  =  -<-Me-Mo> 

r  -  m<Mr-Me-M0>  0<m<  1 

0  "  (m-1) 

r  =  <Mr-M0"MO> 

6  "  -<Mr-M0-MO> 
r  =  g<Mr-M0-Mo)+(l-g)<Mr-Mo> 


A/c0  9  M0  M0^  ^  C  9  ^  1 
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Stress  Complementary 
Regime  Function 


AQB 


Flow  Rule 


A  complete  solution  of  this  plate  problem  consists  of 

a  stress  solution,  M  and  MQ ,  satisfying  the  equilibrium 

IT  y 

equation  (3.7)  and  the  stress  boundary  conditions,  and  a 
velocity  solution  w(r)  which  satisfies  the  velocity  boundary 
conditions  and  is  derivable  from  the  stress  solution  by  use 
of  equations  (3.1)  and  the  flow  rule  of  Table  II. 

The  boundary  conditions  for  this  problem  are 


dw 

dr 

(a)  =0 

(3.10) 

dw 

dr 

(0)  =  0 

(3.11) 

w 

(a)  =  0 

(3.12) 

These  differ  from  those  for  the  corresponding  Tresca 
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In  the  latter 


rigid-perf ectly  plastic  plate  problem 
problem  it  is  possible  to  admit  a  non-zero  slope  at  the 
plate  edge  due  to  the  formation  of  a  hinge  circle.  A  hinge 
circle  is  a  line  along  which  /cQ  is  discontinuous.  Along  a 
hinge  circle  | |  equals  MQ .  If  a  sufficient  number  of  hinge 


32 


circles  form  it  is  kinematically  possible  for  all  the 

bending  to  occur  at  the  hinge  circles  only  and  for  the 

* 

plate  to  collapse  »  This  is  analagous  to  the  formation 

of  plastic  hinges  in  a  beam  or  frame 0  In  the  visco-plastic 

dw 

plate  problem,  discontinuities  in  and  thus  /Cq  are  not 

possible  since  due  to  the  viscous  behavior  of  the  material, 

infinite  stresses  would  be  required  to  produce  such  discon- 
.  .  .  dw 

tinuities.  Consequently  is  continuous  throughout  the 
visco-plastic  plate. 

The  boundary  conditions,  the  equilibrium  equation,  and 
the  flow  rule  constitute  the  complete  boundary  value  problem. 
If  the  entire  plate  lies  in  one  stress  regime  then  the  prob¬ 
lem  is  relatively  simple.  If  however  the  stress  solution 
covers  more  than  one  stress  region  then  each  region  with 
its  corresponding  flow  rule  must  be  considered  separately. 

To  begin  the  solution  an  a  priori  assumption  is  made 
as  to  the  stress  profile.  For  this  problem  it  is  clear 
that  the  stress  profile  includes  more  than  one  stress  regime 
since  must  change  sign  at  some  point  in  the  plate. 

From  circular  symmetry  and  must  be  equal  at  the 
centre  of  the  plate  and  since  compression  occurs  in  the  top 
surface  of  the  plate  at  r  =  0  when  it  deflects  downward, 


*  Plastic  collapse  of  a  rigid-perf ectly  plastic  structure 
occurs  at  that  point  in  the  loading  program  at  which  de¬ 
formation  first  becomes  possible.  The  corresponding  load 
is  called  the  limit  load. 
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both  M  and  M  are  positive  at  that  point.  It  is  reasonable 

¥  xy 

to  assume  therefore  that  in  some  finite  region  around  r  =  0 
the  stress  points  would  lie  in  stress  regime  OB.  Knowing 
also  that  must  be  negative  at  r  =  a,  the  stress  profile 
in  Table  III  is  assumed. 


TABLE  III 

ASSUMED  STRESS  PROFILE  FOR  VISCO-PLASTIC  PLATE  PROBLEM 


Plate 

Stress 

Region 

Regime 

0 

<  r  <  e 

OB 

e 

<  r  <d 

BOC 

d 

<  r  <b 

COD 

b 

<  r  ^  a 

OD 

The  flow  rule  for  each  of  these  regions  in  the  plate 
are  given  in  Table  IV. 

TABLE  IV 

FLOW  RULE  FOR  VISCO-PLASTIC  PLATE  PROBLEM 


Plate 

Region 


0  <  r  <  e 

A/cr 

e  <  r  <  d 

A/cr 

d  <  r  <.  b 

A>cr 

b  <  r  ^  a 

A/Cr 

Curvature  Rates 

f(Mr-M0^  A/Cq  =  (l-f)<Me-M0^>  0<  f  <0.5 

f=0.5atr=0 
f  =  0  at  r  =  e 

0  A/cg  =  <Mg— M0> 

-<Me-Mr-Mo>  A<ce  -  <Me-Mr-M0> 

-t<-Mr-M0)+(t-l)  (Me-Mr-M^  Akq  =  (l-t)(Me-Mr-M0) 

0  <  t  <  1 
t  =  1  at  r  =  a 


■ 


vi  aiaAT 
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Each  of  these  four  regions  in  the  plate  are  now  con¬ 
sidered  separately. 

In  the  first  region,  0  <  r  the  two  bending  moments 

are  equal  and  positive.  Addition  of  the  curvature  rates  gives 


M/Cr+/C9^ 


-A 


Ad/  dw\ 
r  dr  y  dr/ 


n(Mr-M0)  +  (l-n)  . 


Since  M 

r 


A 

r 


the  above  becomes 


(3.13) 


The  equilibrium  equation  reduces  to 


dM 
_ _ r 

dr 


_  ££ 
2 


Integration  gives 
M 


-  + 


=  M 


© 


0  <  r  <  e 


Substitution  of  this  expression  for  M  into  equation 


(3.13)  gives 

A  _d_ 
r  dr 


*  f )  -  -  *§-  +  ci  -  Mo  • 


Integrating  this  result  and  making  use  of  boundary 
condition  (3.11)  gives  an  expression  for  the  velocity. 


4  M,_r‘ 

Aw  -  *§4  +  ^ 


clr 


4-  c. 


0  <C  r  <  e 


* 


The  subscripted  c " s  are  all  constants  of  integration. 
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In  the  second  region  of  the  plate,  e  <l  r  ^d,  the  curva- 

'  d2w 

ture  rate  k  is  seen  to  be  zero.  Therefore  — ~  =  0.  Xnte- 

dr 

grating  this  twice  results  in 


w  =  c3r  +  c4 


e  C  r  si  d 


Using  this  result  and  the  relation  (3,1)  yields 

Ac. 


A/c 


0 


M  -  M 
0  0 


A  dw 
r  dr 


Consequently 


Ac. 


M 


0 


+  M 


0 


e  <1  r  d 


Substitution  of  this  result  into  the  equilibrium  equation 
and  integrating  yields 


M  =  -  +  M 

r  6  0 


XC3  .  .  C5 

■ — —  In  r  +  • — 
r  r 


<  r  <:  d 


In  the  third  region  of  the  plate,  d  ^  r  ^b,  the  curva¬ 
ture  rates  are  equal  in  magnitude  and  opposite  in  sign,  con¬ 


sequently 


n2 

d  w 


dr 


_  .  1.  dw 

2  r  dr 


=  0 


This  is  an  Euler  type  differential  equation  and  is  solved 
by  making  the  substitution  r  =  efc .  The  solution  is 


w  =  o.r  In  r  +  c_ 
6  / 


d 


Substituting  this  into  the  flow  rule  for  this  region  gives 

Ac,. 


_  _A  dw 
*0  r  dr 


—  =  m  -  M  -  M 

2  U©  r  0 


( 3  „ 


Solving  this  for  and  substituting  into  the  equilibrium 


14) 


■  i  i  <  •  • 
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equation  yields 
dM 


~~  +  M 
dr  r 


Ac, 


■^  +  M  +  M„ 
2  r  0 


+ 


£2L 


=  0 


Integrating  this  result  gives  an  expression  for  M 
from  which  Mg  can  be  determined  using  (3.14) .  Thus 

2  Acfi 

Mr  =  M0  In  r  _££-  +  --§+  c8 

2r 


d  ^  ,r  ttC,  b 


M„ 


MQ (1  +  In  r) 


SL. 

4 


Ac, 


2r 


2  +  C8 


In  the  outermost  region  of  the  plate ,  bs^r^a5  the 
bending  moment  M  is  zero.  From  this  and  the  flow  rule 


it  follows  that 


Afc 


-,2 

d  w 

,  2 
dr 


M  +  Nh 
r  0 


Setting  Mg  to  zero  in  the  equilibrium  equation  and 


integrating  results  in 

M  =  - 


2 

£I_ 

6 


b  <r  <a 


Combining  this  with  equation  (3.15)  and  integrating 


gives 


\  _  M  r-  -j-  £r- —  Q  1  ri  r  +  c 

A  dr  w0  18  Cg  in  r  -r  c1Q 


Boundary  condition  (3.10)  determines  the  value  of  c 

3 

c10  =  M0a  “  P18  c9  ln  a 


10 


Integration  of  (3.16)  and  use  of  boundary  condition  (3.12) 


gives 


M, 


Aw  =  -  (r-a)  2  +  (r^-4a^r+3a4)  +  cg(r  ln  —  +  a  -r) 


3.15) 


3.16) 


37 


For  each  of  the  four  regions  in  the  plate  an  expression 
has  been  found  for  the  two  moments  M  and  IVL  ,  and  for  the 
velocity  w.  It  remains  now  to  eliminate  the  nine  constants 
of  integration. 

From  equilibrium  considerations  M  must  be  a  continuous 

function  of  the  radius.  Similarly  Mg,  w,  and  are  also 

continuous  functions  of  r.  Therefore  at  each  of  the  three 

boundaries  separating  the  four  regions  of  the  plate,  these 

continuity  conditions  provide  four  equations. 

dw 

At  r  =  e  continuity  of  w,  and  Mg  yields  the 

following  equations  respectively. 


2 

64 


3 

16 


2 

4 


+ 


2  2 
c  e  Me 

4~~  +  ~4  +  c2 


Xc3e  +  Xc4 


flf  +  V 


=  Ac. 


Ac. 


—  In  e  +  M. 
e  0 


2 

££_  +  _5 
6  e 


(3.17) 


(3.19) 


+  M 


0 


(3.20) 


At  the  boundary  r  =  d  the  continuity  conditions  result 
in  the  following  four  equations. 


c-.d  +  c„  =  cr  In  d  +  c_, 
3  4  6  7 


(3.21) 


c6/d 


3.22) 


Ac. 


-r—  In  d  +  Mn 
d  0 


-.2  c. 

+  __5 

6  d 


Ac, 


2d 


5  +  Mq  in  d 


,2 

P.d._f  c 

4  ^  8 


(3.23: 


i  yssknuod  9ff: 
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Ac. 


Ac 


— —  +  M 
d  0 


2d 


|  +  Mq(1  +  In  d)  -  +  Cg 


The  continuity  conditions  at  r  =  b  result  in  another 
four  equations. 


-M, 


ACg  In  b  +  Ac^  =  (b-a)  3  +  -^“-(b 


-E-fh4  - 


3  4 

4a  b  +  3a  ) 


(3.24) 


+  c_  (b  In  —  +  a  -  b) 
9  a 


(3.25) 


Ac. 


-  Mq  (b-a)  +  -^(b3  -  a3)  +  c9  In  ^ 


(3.26) 


Ac, 


2b“ 


+  Mq  In  b 


-Ac, 


_ 


+  c 


8 


_  £bf  _  fi 

6  b 


0  =  - ^  +  Mq(1  +  In  b)  -  - f  c 


2b 


8 


(3.27) 


(3.28) 


Thus  there  are  a  total  of  twelve  independent  simultaneous 
algebraic  equations  from  which  the  twelve  unknowns  can  be 
found.  The  unknowns  are  the  nine  constants  of  integration 
and  the  three  boundary  radii  e,  d,  and  b. 

The  nine  constants  of  integration  are  first  solved 
for  in  terms  of  the  other  three  unknowns. 

Solving  equation  (3.20)  for  c^  and  substituting  the 
result  into  equation  (3.18)  gives 

c  =  -  2^- 

3  8  A 


 3pe~ 


8 


+  M 


0 


and 


■ 
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From  equation  (3.22) 


3  , 
£e_d 

8  A 


Solving  equation  (3.19)  for  c..  and  substitution  of 

o 

c^  and  gives 

3 

c5  =  £§5  (7  -  3  In  e)  . 

Substitution  of  into  equation  (3.28)  results  in 


'8 


.E. 


16b 


2  (4b4  -  e3d)  -  MQ  (1  +  In  d) 


The  constant  c^  may  be  found  in  terms  of  and 
from  equation  (3.27)  . 


24b  (3e3d  -  4b4)  +  M0b 


Solving  equation  (3.25)  for  results  in 


3  d  In  b  ,  b 

G  * - -  T*  - - 

8  72 


a3b  ,  a4  ,  /e3d  b4\  /  b  a  ' 
T5  +  24  +  (“8"  -  — )  (ln  I  +  b  -  1 


Substitution  of  ,  c^,  and  into  equation  (3.21) 
yields 


M0  /  a2  ,  ,  ,  3i2  ,  ,2  ,  b\  ,  p 

c4  ~  (  ~  ~2~  4  2ab  ~  2~  +  b  ln  I  )  +  A 


13bl  e^d  £ 

72  8  d 


3,  4  /  3  , 

a  h  +  a_  +  fe  a 


18  24 


8 


V)  (  ln  a  +  t 


Finally  from  equation  (3.17)  an  expression  can  be 


2* 


found  for  c 


40 


c 


2 


2 

y  +  2ab 


a3b  .  e3d 

“T8  +  “8" 


+  P 


4  4 

3e  ,  13b 

64  72 


All  nine  constants  of  integration  have  thus  been  found 
in  terms  of  the  radii  e,  d,  and  b.  There  remain  three 
equations,  (3.23),  (3.24),  and  (3.26),  which  have  not  been 

used.  Substitution  of  the  expressions  for  the  constants 
c^  to  Cg  into  these  three  equations  results  in  three  trans- 
endental  equations  from  which  e,  d,  and  b  may  be  determined. 

The  following  non-dimensional  quantities  are  now 
introduced . 


Substitution  for  c^ ,  c^,  c^,  and  c^  into  equation 
(3.23)  and  use  of  the  above  non-dimensional  parameters 
results  in  the  first  of  the  three  equations. 


4 


+ 


e3  5 
16(3  ‘ 


0 


(3.29) 


The  other  two  equations  which  are  found  from  equations 
(3.24)  and  (3.26)  are 


16A  (l  -  In  |)  +  ^  +  462  ~  4P2  +  - 


0 


(3.30) 


and 


(i  +  in  p)  -  fg  (i  -  3  in  p)  +  ^  +  A 


In  (3 


•1  +  p 


0 


(3.31) 
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M 

0 

The  quantity  A  =  „  will  be  called  the  loading  para- 

pa 

meter  since  it  is  a  measure  of  the  loading  on  the  plate 
relative  to  the  yield  strength.  Hopkins  and  Prager 
have  shown  that  for  a  Tresca  rigid-perf ectly  plastic  clamped 
circular  plate  loaded  with  a  uniform  pressure  p,  deformation 
does  not  occur  until  the  pressure  p  reaches  the  critical 


11.26  Mq 

value  pcr  = — ____ _  ^  then  the  plate  collapses.  The  rigid™ 

a 

perfectly  plastic  plate  cannot  support  a  pressure  greater 
than  this  critical  value  without  undergoing  finite  deforma¬ 
tion  which  would  result  in  membrane  stresses.  A  visco™ 
plastic  plate  however  can  support  loads  greater  than  p 

C  JC 

due  to  viscous  stresses.  Therefore  in  the  visco-plastic 
plate  problem  the  loading  parameter  A  can  take  values  from 
1/11.26,  corresponding  to  the  special  case  of  a  rigid-perf ectly 
plastic  plate,  to  the  value  zero  representing  a  visco-plastic 
plate  with  zero  yield  moment. 

Hopkins  and  Prager  found  also  that  for  the  Tresca  rigid- 
perfectly  plastic  plate  the  boundaries  between  the  stress 
regions  are 


e  =  0,  6  =  0.730,  p  -  1.0 


Substitution  of  these  values  of  e,  6,  and  f3  into  equa¬ 
tions  (3.29),  (3.30),  and  (3.31)  along  with  A  “  1/11.26  shows 

that  this  solution  thus  far  reduces  to  that  found  by  Hopkins 
and  Prager  for  the  special  case  of  the  rigid-perf ectly  plastic 
plate . 
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Solutions  for  €,  6,  and  f3  have  been  found  for  a  number 
of  values  of  A  ranging  from  1/11 „ 26  to  zero.  These  results 
have  been  found  with  the  aid  of  a  computer  using  an  itera- 
tive  method  based  on  Newton's  method  for  obtaining  the  roots 


of  an  equation 


11 


The  variation  of  the  three  stress  region 
boundaries  with  the  parameter  A  is  shown  in  Fig „  3 . 3 „ 

By  substituting  the  solutions  for  the  nine  constants 
c-^  to  Cg  into  the  equations  for  w ,  M  ,  and  Mg  for  each  of 
the  four  stress  regions,  a  complete  solution  to  this  problem 
is  obtained,  since  for  any  value  of  the  loading  parameter 
A  the  stress  and  velocity  solutions  are  completely  determined 
with  the  aid  of  Fig,  3.3. 

The  stress  solution  thus  obtained  is 


M 

M 


r  _ 


=  1  + 


0 


M 

M 


M 

i 

M 


=  1  + 


=  1  + 


8A 


8A 


,  2  0  2 
3e  -  2p 


_2  0  2 
3e  -  2 p 


0 


M 

M 


1  +  i 


7e  p  e  In  e/p 
_24p  6  8p 


€3/8p 


M 

M 


r  _ 
0 


In  -  1  +  — 
P  A 


1,Q2  2\ 
-  p  ) 


M 

i 

M 


£  = 
0 


ln  f  +  i 


'S 


0  ^  p  £ 


f  €  <  p  <  6 


els  n  +  _i_ 

16  V  2  ~  ft2 


1  ,ft2  2*  ,  €  6  / 1 

4  (P  -  P  )  +  ~I6  (  — 

VP 


r  6 
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M 
_ 2 

M 


0 


-£  + A 

p  A 


2  3  r  q3‘ 

_  £ _ _§ 6  ,  jr 

6  8Pp  +  6p 


M 

M 


®  =  o 
0 


P  <  p  <  1.0 


The  complete  velocity  solution  is 


wA 

Moa' 


+  2P 


3 £ 


+  p"  In  p 


•4  .22 

£ -6J..  .P. 


-  3e 


64 


+ 


+  ln 


72 


8 


£  .  £  +  _i_  +  (i±  .  Si 

5  18  24  l  8  6 


1"  P  +  p 


0  <  p  <C  e 


wA 

M0a: 


1  +  2p  -  ^ 

2  2 


+  p  In  P 


3 

: _ i 

8 


-  Al£  +  ^  In  £ 


8 


+  12&L  _  JL  +  -L  + 

72  18  24 


if*  j£ 

8  6 


ln  P  + 


P 


p  ^  6 


wA  _ 
2 


Moa 


1  +  2P  -  ^ 

2  2 


+  p^  ln  p 


3  4 

O  m  £  +  £ 

8  n  p  72 


AL  -p  JL - p 

18  24 


3  <•  «4\  ✓ 


£6 

8  "  6 


Aln  p  +  t 


-  1 


6  ^  p  ^  p 


wA  


Moa 


P(1  -  p)  +  Pp  ln  p  -  J  (p  -  1) 


!]*J 


r  4 
£_ 
72 


+ 


•As  £' 

8p  “  6 


p  In  p  +  1  -  p 


18  24 


P  <  p  1.0 


When  the  viscous  stresses  in  the  plate  become  very 


small  in  comparison  with  the  stresses  causing  plastic 
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deformation  this  solution  approaches  that  found  by  Hopkins 
and  Prager  9  for  the  corresponding  Tresca  rigid-perf ectly 
plastic  plate. 

It  may  be  verified  that  parameters  f  and  t  as  defined 
in  Table  IV  lie  in  the  following  closed  intervals. 


0  <  f  <0.5 


and 


0 


Also  it  may  be  shown  that 


f  =  0 , 5  at  r  =  0 


t  =  0  at 


b 


f  =  0 


at 


t  =  1  at 


Numerical  results  for  the  stress  and  velocity  solution 
have  been  obtained  for  values  of  A  ranging  from  1/11.26  to 
zero  and  are  plotted  in  Figures  (3.3),  (3.4),  (3.5),  (3.6), 

(3.7)  and  (3.8) „ 
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FIG.  3.5 

VISCO-PLASTIC  PLATE 


STRESS  REGION  BOUNDARIES  vs.  LOADING  PARAMETER 
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FIG.  3.4 

VISCO-PLASTIC  PLATE 
DEFLECTION  RATE  vs.  RADIUS 


I- 


t 


' 

-  *.•* 
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1.0 


FIG.  3.5 

VISCO-PLASTIC  PLATE 

B  EM  DING  MOMENT  M  vs.  RADIUS 

r 
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FIG.  3.6 

VISCO-PLASTIC  PLATE 
BENDING  MOMENT  M„  vs.  RADIUS 
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FIG.  3.7 
VISCOUS  PLATE 


DEFLECTION  RATE  vs.  RADIUS 


. 


d^/W  1N3W0W 
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CM 


FIG-,  3.8 
VISCOUS  FLATE 


BINDING-  MOMENT  vs.  RADIUS 


CHAPTER  IV 


CREEP  PROBLEM 


A  different  constitutive  relation  is  now  applied  to 
a  sandwich  plate  otherwise  identical  with  that  in  Chapter  III 
The  method  of  solution  which  is  used  was  first  introduced 


by  Venkatraman  and  Hodge 


16 


however  the  complementary 


function  was  not  introduced.  The  method  was  applied  to 
both  simply  supported  and  clamped  circular  plates. 

The  constitutive  equation  for  the  material  now  consid¬ 
ered  can  be  obtained  from  the  complementary  function  which 
is  assumed  equal  to  a  function  of  the  generalized  stress  in¬ 
tensity  M  where 


M  =  max 


M 


M, 


M 


M. 


0 


(4 


times  a  function  of  time.  If  stationary  creep  only  is  con- 

(r  ■' 

sidered  then  Ec  may  be  chosen  as 


E 

c 


(M)  n+1 

(n+1)  |j,n 


G  ( t ) 


where  and  n  are  creep  constants  depending  on  the  material 
and  the  temperature  and  jx  has  the  dimensions  of  moment  per 
unit  length.  The  constitutive  relation  found  from  (4,2) 
represents  time  hardening  creep.  The  time  hardening  effect 
is  introduced  by  the  term  G(t) .  This  surface  of  constant  Ec 
is  a  hexagonal  prism  in  principal  stress  space. 
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The  principal  curvature  rates  are 

bE  dE 

Kr  =  and  Ke  =  cTmT  -  (4‘ 

r  9 

When  n  approaches  infinity  the  constitutive  relation 
defined  by  (4.2)  and  (4.3)  approaches  that  for  the  Tresca 
rigid-perf ectly  plastic  plate  problem  with  a  yield  moment 
equal  to  |_i.  When  n  is  one  this  relation  is  essentially  equiva¬ 
lent  to  that  given  by  (1.23). 

The  complete  flow  rule  derived  from  (4.2)  using  (4.3) 
and  the  condition  that  the  rate  of  energy  dissipation  per 
unit  area  is  constant  on  the  Ec  surface  is  given  in  Table  V. 


TABLE  V 

FLOW  RULE  FROM  EQUATION  (4.2) 


Stress 

Regime 


COB 


OC 


Complementary 

Function 


M 


n+1 


9 


(n+1)  [i 


n 


G(t) 


M 


n+1 


0 


(n+1)  \i 


n 


G  (t ) 


Flow  Rule 


jc  =  0 


K 


9 


K.  ~ 


K-rs  ~ 


'V)n 

v  M-y 


G  ( t ) 


"M0"Mr\n 

-s  I  — - -  )  G  ( t ) 


0  ^  s  1 


+U\  n 


£Y 


'9 


(1-s)  f  —  j  G(t)  +  s 


/M.-M  \n 

,-v) 


G  ( t ) 
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Stress  Complementary 
Regime  Function 


Flow  Rule 


COD 


(Me-Mr)rri'1G(t) 


(n+1)  n 


n 


Kr  =  -(-V)  G(t) 


/Mq-M  \n 

'co-  +  (-V)  G(t) 


OD 


(-M  ) 
r 


n+1 


(n+1) 


n 


G(t) 


/-M  \n  /Mo“M  \n 

Kr  =  (-/)  G(t)  +  — JT^1)  G(t) 


/Mg  M  \n 

/Cg  =  (l-t)f jjj - — J  G (t ) 


0  <  t  <  1 


DOE 


( -M  ) 
r' 


n+1 


(n+1)  ii 


n 


G(t) 


/-M  \n 

Kr  ■  -m  ^ 


K9  =  0 


OE 


(~Mg) 


n+1 


(n+1)  (i 


n 


G  ( t ) 


-  -  (+)"-> 


Ke  =  (<j-1)(it)  G(t) 


0  <  q  <  1 


EOF 


n+1 


(n+1)  |i 


n 


G  ( t) 


/c  =  0 
r 


K 


0 


.-Mgxn 

.  U 


G  ( t) 


OF 


(>Mg) 


n+1 


(n+1)  |i 


n 


G  ( t ) 


/M  -Mg\n 

*r  "  <^)  G<t> 


/cg  =  ( m- 1 )  G(t) 


0  <  m  <r  1 


/M^”Mg\n 


G(t) 
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Stress 

Regime 


FOA 


OA 


AB 


Complementary 

Function 


(Mr-M9) 


n+1 


(n+1)  \i 


n 


G(t) 


(M  ) 
r 


n+1 


(n+1)  [i 


n 


G(t) 


(M 


n+1 


(n+1)  n 


n 


G  ( t ) 


Flow  Rule 


K 


/Mr”M0xn 


G  (t ) 


,Mr-Mexn 

G(t) 


M  -M  \  n  M  \  n 

Kr  -  g(-Vv G(t)  +  (1-g)(lf)  G(t) 


Ke  =  G(t) 


u 


/M  \n 

Kr  -  uf) G(t) 


0  g  <C  1 


OB 


(M  ) 
r 


n+1 


(n+1)  [i 


n 


G  (t ) 


*9  =  ° 


K 


/M\n 

=f(ir 


G  (t ) 


1C 


@ 


0  <  f  <  1 


An  a  priori  assumption  is  made  for  the  stress  profile 
The  stress  profile  of  Table  III  is  assumed. 

The  boundary  conditions  for  the  problem  are 


dw 
dr 

dw 
dr 

w(a) 


(0) 

(a) 


=  0 

=  0 
=  0 


(4.4) 

(4.5) 

(4.6) 
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As  in  the  previous  problem  the  bending  moments  must 
satisfy  the  equilibrium  equation  (3,5)  and  the  curvature 
rates  must  be  derivable  from  the  velocity  using  equations 
(3.1)  . 


M 


In  region  OB  the  stress  intensity  M  is  equal  to  and 
g.  Therefore  addition  of  the  curvature  rates  gives 

/Me\n 

/cr  +  Kg  =  (~J  G(t)  .  0<  r  O 


(4.7; 


The  equilibrium  equation  in  this  region,  0  r  e , 


becomes 


dM 


r  +  EJL 


dr 


0 


Integration  of  this  gives 


m  =  —  ££ —  -j-  c  —  m 

r  4  C1  ©  * 


(4.8) 


From  equations  (3.1),  (4.7),  and  (4.8) 


1  d  f  =  (-^}  G  (t )  =  G  ( t ) 


r  dr 


dr 


n 


Integration  and  the  boundary  condition  (4.4)  gives 


dw 

dr 


2G 


PHn(n+l) 


2JL 


+  Cl) 


\  n+1 


n+1 


from  which 


w  = 


2G 


PUin(n+l) 


2  \  n+1 

+  c  ) 

4  C1 J 


-  c. 


n+1 


dr 

r 


(4.9! 
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In  the  second  region  of  the  plate,  e  ^  r  ^  d,  the  curva¬ 


ture  rate  ic  is  zero  and  therefore 


w  =  -c2  r  +  c3 


e  <  r  C  d 


(4,10) 


From  Table  V  and  equation  (3.1) 


which  gives 


M 


0 


11 

r 


'C2\1/n 


Gr 


(4.11) 


Substituting  this  into  the  equilibrium  equation  and 
integrating  results  in 


M 


_m_  ft i\1/n  _  esI  +  t± 

(n-1)  VGry  6  r 


e  ^  r  ^  d 


This  expression  is  valid  for  n  not  equal  to  one 
When  n  is  equal  to  one 


M 


Lie  _  2c. 

2  pr  ,  4 

-77-  In  r  -  •£~ - r  — 

Gr  6  r 


This  latter  expression  will  not  be  used  thereby  restrict 
ing  n  to  be  other  than  one.  The  particular  problem  when  n 
does  equal  one  is  equivalent  to  the  visco-plastic  plate  prob¬ 
lem  solved  in  Chapter  III,  with  the  yield  moment  equal  to 
zero . 

In  the  third  region  of  the  plate,  d  ^  r  ^b,  the  curva¬ 
ture  rates  are  equal  in  magnitude  and  opposite  in  sign  so 
that  addition  gives 
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a2 
d  w 

,  2 
dr 


dw 

dr 


=  0 


from  which 


w 


•c.  In  r  +  cr 
5  6 


d  <  r  <  b 


Use  of  (4.13)  in  the  egression  for  k.^  in  this  region 


from  Table  V  gives 


or 


_5 

2 


M.  -  M 
0  r 


4VL  -  M  \n 
©  r  \ 


|l 


M* 


C5  \l/n 


-Gr 


Substitution  of  this  latter  expression  into  the 
equilibrium  equation  gives 


dM 


dr 


M- 


\l/n  2 

2 )  2 


=  0 


from  which 


M 


-m 

2 


-Gr 


'Cr  \l/n 


-Gr' 


< 


££-  +  c 
4  7 


1 


and 


M@  =  u 


2-n 


1/n 


-Gr‘ 


££__  +  c 
4  7 


►  d 


(4  0 12) 


4.13' 


(4.14) 


(4.15) 


(4.16) 


(4.17) 


In  the  outermost  region  of  the  plate,  b  <C  r  _<c  a,  the 
moment  is  zero.  Thus  integration  of  the  equilibrium 
equation  gives 


M 


2  cQ 

+  __8 

6  r 


(4.18) 


From  equation  (4.18)  and  the  flow  rule  for  this  region 


d2w 

dr2 


/  M  \n  /  2  c._\n 

-  -<-t)  ■ 


(4.19) 


:  .  M  Jnsoiom 


58 


Integration  of  this  and  use  of  boundary  condition  (4.5) 


gives 


dw 

dr 


d 


n 


%a  r  „2 

PC  c 

6 


8 


c 


n 


(4.20) 


Integrating  once  more  and  use  of  boundary  condition  (4.6) 
results  in 
G 


w  = 


,a  ^  ^2 


c8\n 


n  ,  i  I* 


c 


dCd| 


(4.21) 


Equilibrium  requires  that  and  Mg  be  continuous.  Also 

.  dw 

the  viscosity  effect  requires  that  be  continuous,  and  since 

the  plate  does  not  fracture,  w  is  continuous.  Therefore  these 

continuity  conditions  applied  to  the  three  boundaries  of  the 

stress  regions  result  in  a  total  of  thirteen  equations. 

dw 


Continuity  of  w,  M^,  and  Mg  respectively  at  r  =  e 


gives 


2G 


n 


PH  (n+1) 


-  ^-+  Cl 


.n+1 


-c. 


n+1 


dr 

r 


i 


r=e 


-  -c^e+c^  (4.22) 


2G 


PHn (n+1) 


 £§. 


n+1 


+  c. 


-  c. 


n+1 


(4.23) 


J2§_  +  c 

4  C1 


=  BL 


n-1  \Ge 


2^\1/n  .  £e!  +  f4 


(4.24) 


 22. 


+  c. 


d 


^y/n 

GeJ 


(4.25) 


At  r  =  d  these  continuity  conditions  and  the  condition 
Mr(d)  =0  gives 


•n 
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-  c^d  +  =  -Cj.  In  d  + 


-  c. 


-c5/d 


m-  (Hi Y/n  .  m!  +  fi  , 

n-1  \  Gdy  6  d  u 


_  m  (c  s  -  M? 

2  \  „  -.2  /  4 


'•Gd 


+  =  0 


M- 


c2\1/n 


Gdy 


n 


c  ^  \l/n 

V 


’2-nVv5  V'  “  _  jod 


-Gd 


4  +  °7 


At  r  =  b  the  following  equations  are  obtained 


(4.26) 

(4.27) 

(4.28) 

(4.29) 

(4.30) 


-Cj.  In  b  +  c,  =  — 
5  6  n 

d 


-c,/b  =  - 
5  n 


P  a  p  a 

-  *2  c0-i 

PC  8 

4  \ 

L 6  e  _ 

dCd^ 


nu  /c5  \i/n  ~  pb2 
2  V  ^  .2  )  4 


y 

r  r2  C0n 

Ei_  8 

L 

L 6  t, . 

n 


dC 


+  c. 


-Gb 


,2  c0 

+  _8 

6  b 


*  (•’?)  ©v"  -  *  * • • 


(4.31) 


(4.32) 


(4.33) 


(4.34) 


The  unknowns  are  e,  d,  b,  to  Cg,  and  one  other  constant 
of  integration  from  equation  (4.9).  There  is  a  total  of  twelve 
unknowns  and  what  appears  to  be  thirteen  simultaneous  independ¬ 
ent  equations.  All  thirteen  equations  are  not  independent 
however  since  equation  (4.29)  is  a  combination  of  equations 
(4.27)  and  (4.30) . 

Using  equation  (4.22)  to  eliminate  the  constant  of 
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integration  in  equation  (4.9)  gives 


w 


2G 


n 


n+1 


P|U  (n+1)  Jr  L 


2  \n+l 

+  c 

4  1 


C 


C2e  +  C3 


0  ^  r  ^  e 


:4.3s: 


Solving  equation  (4.28)  for  gives 

c  =  £d!  .  md  fZ 2\1/n 
4  6  n-1  \GdJ 


4.36) 


in 


Subtracting  equation  (4.25)  from  equation  (4.24)  results 

n-1 

+  2\1/n  [YA1/n  _  nd  n  1 

1  V  <3  )  W 


*  (*3-3)  *  £ 


0 


4.37) 


Solving  equations  (4.36)  and  (4.37)  for  and 
yields 


c2  ~  G 


p  (e3-d3)  (n-1)  1  n 
n-1  n-1 


6|i^e 


-  nd 


and 


n-1 


=  2S L  _ 


p(e3-d3) (n-1) 

n-1  n-1 


6  e 


n 


nd 


n 


Substituting  c_  in  equation  (4.25)  results  in 


2  /  3  ,3W  -i  x 

r  _ _ p (e  -d  ) (n-1)  __ 

1  4  n-1  n-1 

6e1/n  fe  n  -  nd  n 
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From  equation  (4.27) 


Gd 


p (e3-d3) (n-1) 

n-1  n-1 


n 


6  |-l  e 


n 


nd 


n 


The  above  and  equation  (4.29)  or  equations  (4.27)  and 
(4.30)  gives 

_  n 


12 


p(e3-d3) (n-1) 


L 


n-1  n-1 

d1/n  (e  n  -  nd  n 


2 


+ 


4  * 


Substitution  of  the  expressions  for  c,.  and  into 
equation  (4.33)  yields 

n—  2 

c  =  J  -  IsL  +  +  njnzjlE(e3-d3.L-  (~±-  _  b”  dVrA 

8  P\  12  4  )  n^l  n-1  1/n  b  d  / 

12 ( e  n  -  nd~n“^ 


From  equation  (4.31)  the  following  is  obtained 


c,  =  Gd  In  b 
6 


p (e3~d3) (n-1) 

n-1  n-1 


n 


^a  pa 


+ 


6|i  (  e  n  -  nd  n 


d 


n 


i 


-E-  »  2^2  +  b 


12 


c 


n-2 


3bd2\  _  n  (n-1)  p  (e3-d3)  (b/dly/n  -  b  n  d1/n)~|n 


n-1  n-1 

12  f  e  n  -  nd  n 


dCdi 


Equation  (4.26)  gives 


Gd(l+ln  -g) 


p(e3-d3) (n-1) 

n-1  n-1 

6| i  (e  n  -  nd  n 


n 


n  ,  /-u  V  It 

d  Jb 


a  r,  a  r  /  ~  ,  3 

f 9^2  b_ 
12  1  ^  C 


n-2 

3bd2\  n  (n-1)  p  (e3-d3)  (b/d1/n  -  b  n  d1/n_in 


C 


n-1  n-1 

12  ( e  n  -  nd  n 


d^d|  . 
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The  three  remaining  equations  may  now  be  used  to  find 
the  values  of  e,  d,  and  b„  The  following  non-dimensional 
parameters  are  used* 


Substitution  of  the  expressions  obtained  for  the  constants 
c ^  to  Cq  into  equations  (4.23),  (4.32),  and  (4.34)  gives  the 

following  three  equations  respectively, 


(n+1) € 


(e3-53) (n-1) 


n-1 

6€1/n  U  "  ■ 


n6 


n-1 

n 


n+1 


(€3-53) (n-1) 


6e 


1/n 


n-1 

n 


1  n+1' 
n-1  1 


■n5 


n 


j 


(e3-63)  t n-1) 


n-1  n-1 


(4.38) 


2(e3-S3) (n-1) 


n-1  n-1 


2C2  + 


P3/C  -  3p62/c 


n (n-1)  (e3-63) 


n-1  n-1 


n-2 


=  0 


( 4  0  39) 


(e3-63) (n-1) 


n-1  n-1 


(2-n) 


1/n 


n 


& 


1/n 


+  3(6  -P 


2, 


«  0 


4.40) 
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These  three  Equations  must  be  solved  simultaneously 
for  each  value  of  n.  Solutions  have  been  found  using  a 
computer.  The  integration  in  equation  (4.39)  was  done  num¬ 
erically  using  a  24  point  Gauss  Quadrature  routine  11 
The  roots  of  the  three  equations  were  found  using  Newton's 
Method.  The  results  are  shown  in  Table  VI. 


TABLE  VI 

STRESS  REGION  BOUNDARIES  FOR  CREEP  PROBLEM 


n 

e 

6 

13 

2 

0.43009 

0.67414 

0.80074 

4 

0 . 3!2217 

0.68378 

0.85900 

6 

0.26745 

0.69077 

0.88934 

8 

0.23332 

0.69581 

0.90820 

10 

0.20951 

0.69959 

0.92116 

12 

0.19170 

0.70254 

0 . 93066 

The  values  of  e,  6,  and  (3  for  n  =  1  are  the  same  as 
those  for  the  visco-plastic  problem  when  the  yield  stress 
goes  to  zero. 

As  n  approaches  infinity  the  constitutive  relation 
from  equations  (4.2)  and  (4.3)  approaches  that  for  a  Tresca 
rigid-perf ectly  plastic  plate.  The  results  in  Table  V  show 
that  as  n  becomes  large  the  values  of  6,  and  (3  approach 
rather  slowly  those  for  the  Tresca  rigid-perf ectly  plastic 
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plate  problem*  The  values  for  the  latter  problem  being 
e  =  0,  6  =  0.730,  and  (3  =  1.0. 

Substitution  of  the  values  of  to  into  the  ex¬ 
pressions  for  w,  M^,  and  for  each  region  of  the  plate 
results  in  a  complete  solution  since  e,  6,  and  f3  can  be 
determined  for  any  value  of  n  from  equations  (4.38),  (4.39), 

and  (4.40) . 

The  velocity  solution  is 


w  =  Ga 


2  {2*1 


2\n 


6n 


1 


3e 


+ 


(e3-63)  (n-1) 


-in+1 


n-1  n-1 

1/n  /  n  s  n 

:  e  -  n5 


3  ,  2  w2\ 

.  (e  -£  ) 


(e3-  63) (n-1) 


n-1  n-1 

e1/n  (e  n  -  nd  n 


n+1 


■  *£  + 

c 


6  +  61n  p/6  -  € 


(e3-63) (n-1) 

n-1  n-1 


+  B 


n 


h6 


n 


w  —  Ga3  f  ( 6  +  6 In  P/ 6  —  p) 


(€3-53) (n- 


n-1 

n 


n5 


n-1 


n 


+  B  )  €  ^  p  <  6 


w  = 


(€3-63)(n-l) 


n-1  n-1 


6  <  p  ^  P 


' 
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w  —  Ga. 


o  /  2\  n 

2  ( 1 


1  n  1 


6|l/  nJ 


pj  i 


c2  +  fj 


2C 


n(n-l)  (e3 -  5  3 


2  € 


n- 1 
n 


n-2 

1/n  n  -1/n*  1n 

~  P  o  ) 


n5 


n-1 

n 


1  nlr 


Where  B  = 


P  I 


V2  .  P_  366' 

C  +  2?  '  2T 


n-2 

n(n-l)  (g3-63)  (g/51^11  -  g  n  51/n)  ~>n 


2  € 


n-l 

n 


n6 


n-1 


d£d£  . 


The  stress  solution  is 


M  =  Mn 
r  @ 


/  2n 

£|_ 


3  /  £- 2  2, 
7  (e  -p  , 


+ 


( e  3  —  6  3 )  (n-1) 

n-1  n-1 


:1/n(i  n 


0  ^  p  <^7  € 


n6 


n 


M 


=  SSL 

r  6 


2  r/,3  3 


P 


+  < 


1/n 


i  ©if 


(€3°63)n 
n-1  n-1 

:  n  -  no  n 


2 

Mn  = 

9  6 


(n-1) (c3-a3) 

n-1  n-1 

e  n  -  n5  n 


1/n- 


€  ^  p  ^  6 
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(a2-P2) 


61/nY  ( € 3  —  53 )  (n-X)  V 
2/ n  \  n-1  n-1  ) 

P  n  n 


( 62“p2) 


+ 


(e3-63) (n-1) 


n-1  n-1 


6  <  p  ^  p 


M 


■P 


£ 

2p 


4- 


n-2 
p  n 


P  ^  P  ^1 


M@  =  0 


Numerical  results  have  been  computed  for  w,  M^,  and 
Mq  using  the  above  expressions.  These  values  have  been 
plotted  in  Figures  (4.1),  (4.2)  and  (403)0 
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CREEP  PROBLEM 
DEFLECTION  RATE  vs.  RADIUS 
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FIG.  4.2 
CREEP  PROBLEM 

BENDING  MOMENT  M  vs.  RADIUS 

r 
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CREEP  PROBLEM 

BENDING  MOMENT  M@  vs.  RADIUS 


CHAPTER  V 


CONCLUDING  REMARKS 


Venkatraman  and  Hodge 


16 


have  presented  a  solution  for 


the  stationary  creep  problem  in  Chapter  IV.  That  solution 
is  not  the  same  as  the  solution  found  in  this  thesis.  This 
discrepancy  is  now  discussed. 

Venkatraman  and  Hodge  solved  both  the  simply  supported 
and  built-in  circular  plate  problems  using  the  constitutive 
equation  derivable  from  (4.2).  In  the  former  problem  the 
stress  solution  was  assumed  to  lie  in  stress  regimes  OB  and 
BOC  (referred  to  Fig.  3.2)  whereas  in  the  latter  problem  the 
stress  solution  was  assumed  to  lie  in  regimes  OB  and  OD. 

Consider  the  region  of  the  built-in  plate  whose  outer 
boundary  is  the  circle  along  which  M  is  zero.  The  free 
body  diagram  for  this  portion  of  the  plate  which  is  supported 
entirely  by  vertical  shear  is  identical  with  that  for  a  simply 
supported  plate.  This  requires  that  the  stress  solution  for 
this  inner  region  of  the  clamped  plate  be  the  same  as  that 
for  the  simply  supported  plate  and  the  velocity  solutions 
must  differ  only  by  a  constant.  This  necessary  requirement 
cannot  be  met  using  the  stress  profiles  assumed  by  Venkatraman 
and  Hodge.  Since  a  stress  profile  consisting  of  regions  OB 
and  BOC  was  assumed  for  the  simply  supported  plate,  the  stress 
profile  for  the  clamped  plate  must  also  include  these  two 
stress  regions. 


. 

beJ'ioqque  yJ^' 


71 


If  Venkatraman  and  Hodge  were  able  to  find  a  solution 
for  the  clamped  plate  using  a  stress  profile  consisting  of 
stress  regions  OB  and  OD  then  the  question  arises  can  a 
solution  be  found  for  the  simply  supported  plate  using  a 
stress  profile  consisting  of  stress  region  OB  only?  The 
answer  is  no  and  furthermore  the  solution  found  by  Venkatraman 
and  Hodge  for  the  built-in  plate  can  be  shown  to  be  incorrect. 

If  an  attempt  is  made  to  solve  either  the  visco-plastic 
or  creep  problem  of  a  simply  supported  circular  plate  using 
a  stress  profile  consisting  of  stress  region  OB  only,  expres- 
sions  can  in  fact  be  found  for  w„  M  ,  and  „  Further  in- 
vestigation  however  indicates  that  the  flow  rule  is  violated 
since  it  can  be  shown  that  the  parameter  f  in  Tables  II  and 
V  does  not  lie  in  the  closed  interval  Q  f  ><C  I  but  in  fact 
approaches  infinity  as  r  approaches  a.  Similarly  for  the 
clampled  plate,  if  a  stress  profile  consisting  of  stress 
regions  OB  and  OD  is  assumed,  with  r  ~  a  being  the  boundary 
between  the  two  regions  of  the  plate,  it  can  be  shown  that 
f  approaches  infinity  as  r  approaches  a.  These  solutions 
are  thus  not  admissible, 

Venkatraman  and  Hodge  presented  the  flow  rule  for  stress 
region  OB  as 


M 


r 


M 


and 


£  G(t) 


This  is  equivalent  to  that  shown  in  Table  V,  The  solution 
they  obtained  for  the  clamped  plate  violated  this  flow  rule 
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since  fc  becomes  negative  as  r  approaches  a. 

It  follows  from  this  discussion  and  the  uniqueness  theorem 
in  Chapter  II  that  there  is  only  one  complete  solution  to 
each  plate  problem  using  this  piecewise  linear  form  of  the 
surface  of  constant  Ec .  Any  other  kinematically  admissible 
velocity  solutions  or  statically  admissible  stress  solutions 
which  may  be  found  can  be  shown  by  use  of  the  extremum  prin¬ 
ciples  (2o10)  and  (2„12)  not  to  be  the  real  solutions 0 

There  are  a  number  of  areas  of  research  which  could  follow 
the  work  presented  in  this  thesis „  The  first  of  course  is 
experimental  verification  of  the  two  solutions  using  metals 
near  their  recrystallization  temperatures  so  that  viscous 
effects  are  present,.  Other  problems  include  the  use  of  other 
creep  laws  based  on  the  complementary  function  such  as 
Ec  =  |j,AeM/^  where  |i  and  A  are  creep  constants  depending  on 
the  material  and  temperature  and  M  is  as  defined  in  (401)  „ 

The  problem  of  circular  plates  exhibiting  viscous  effects 
which  are  subjected  to  time  varying  loadings  such  as  a  blast 
or  pulse  loading  is  also  of  technical  importance „  Solutions 
for  this  class  of  problems  would  require  consideration  of 


inertia  effects.. 


V 
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